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ABSTRACT 



Aims. We study the presence of reversed polarity magnetic flux in sunspot penumbra. 

Methods. We applied a new regularized method to deconvolve spectropolarimetric data observed with the spectropolarimeter SP 
onboard Hinode. The new regularization is based on a principal component decomposition of the Stokes profiles. The resulting 
Stokes profiles were inverted to infer the magnetic field vector using SIR. 

Results. We find, for the first time, reversed polarity fields at the border of many bright penumbral filaments in the whole penumbra. 
Key words, methods: observational - methods: numerical - Sun: magnetic topology - sunspots - techniques: polarimetric 



1. Introduction 



MHD simulations of sunspots ( R empel et all l2009allbi : iRe mpel. 
l2012l) have successfully reproduced many aspects of sunspot 
penumbrae as produced by magneto-convection in inclined mag- 
netic fields. As a consequence, magnetic field with inverse po- 
larity should be observed all along the boundaries of bright 
penumbral filaments. T his reversed flux has been observed 
in the outer penumbra dWestendorp" P laza et al 1 [19971 l200ll 
Idel Toro Iniest a~et all 1200 ll) . but even with the highest spatial 
resolution spectropolarimetry , we have been unable to observe 
them in the inner penumbra (Langhansetal. , 2005). When in- 
verting data from the Advanced Sto kes Polarimete r (Elmore, 
1 19921) under the MISMA hypothesis. ISanchez AlmeidaT 72005) 
found that penumbral filaments could be interpreted as bundles 
of magnetic fie ld line s rangi ng over all inclination angles. In 
his PhD thesis, iFranzl d201 ll) shows how the opposite polarity 
fields in penumbral filaments are responsible for a weak third 
lobe in Stokes V that appears on the red wing of 40% of the 
penumbral spectra of the two iron lines at 630 nm observed with 
the Hinode spectrop olarimeter (HINODE/SP, Li tes et all 1200 U 
iKosuei et all l2007h. For rece n t reviews see Borrero & Ichimoto l 
(1201 lh . iBellot Rubiol (1201 Ol) , iBorrerol d2009l) . ISchlichenm aier 
(200$, or lTritschlerfd2009l 



2. Observations 

The active region NOAA 10953 was mapped at an average helio- 
centric angle of 0=12.8° using HINODE/SP on 2007 Apr 30, be- 
tween 18:35 and 22:30 UT. T his active region has recently been 
analyzed bv llchimotol d2010h and lLouis et aU (l2009ll20Tlh . The 
region was scanned in ~1000 steps, with a step width of 07148 
and a slit width of 07158, recording the full Stokes vector of 
the pair of the neutral iron lines at 630 nm with a spectral sam- 
pling of 21.46 mA. The spatial resolution was ~0732. The in- 
tegration time was 4.8 s, resulting in an approximate noise level 
of 1.2 x 10~ 3 . To calibrate the spectra, we averaged the inten- 
sity profiles coming from a 5"x5" region of quiet Sun near the 
disk center. We compared this average profile to the FTS spec- 



tral atlas dKurucz et all Il984t iBrault & Neckel 119871) once it 
was convolved with the spectral point spread function (PSF) of 
Hinode. We have found that a veil of 3.4% must be subtracted 
from the continuum intensity before normalization. This value 
is similar to the stray -light contaminati ons of 4.7% found b y 
iDanilovic etaTI d2008l) and 5% found bv lSocas-Navarrold201 lh . 
The wavelength was calibrated by assuming that the average um- 
bral profile is at rest. 



3. Deconvolution 

Deconvolving two-dimensional spectropolarimetric data wave- 
length by wavelength has some drawbacks. First, the number of 
deconvolutions one has to carry out is large. Second, many of 
these wavelengths contain practically no relevant information in 
Stokes Q, U and V, apart from the noise. This potentially leads 
to an enhancement of the noise level unless it is filtered. For this 
reason, and in order to overcome these two problems, we intro- 
duced a regularization into the deconvolution (the details will be 
given in Asensio Ramos et al. in prep). This regularization acts 
on the spectral dimension of the data by assuming that the orig- 
inal Stokes profiles (before reaching the telescope) at each pixel 
can be written as a truncated linear combination of the principal 
components {(pi(A)} obtained from the observations, so that 



(1) 



where N is the number of eigenfunctions needed to reproduce 
the profiles to adequate precision. Under this assumption, and 
assuming that the monochromatic PSF is independent of wave- 
length, the observed Stokes profiles is given by 



1(A) = £(«,* P) 0,(4) + N, 



(2) 
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Fig. 1. Deconvolved continuum intensity. White arrows corre- 
spond to the horizontal component of the magnetic field in points 
with reversed polarity. Black rectangles mark the areas area am- 
plified in Fig. |3] The black arrow points to the disk center. 

where P is the PSF (kindly provided by M. van NoorQ), and 
N is a Gaussian noise component with zero mean and variance 
cr 2 . Given the orthonormality of the principal components, the 
projection images fulfill: 



<IW),0*W)> = <u**P + N. 



(3) 



where the noise is still Gaussian with zero mean and variance 
cr 2 . Consequently, the regularization process we have used im- 
plies that we have to deconvolve the projections of the origi- 
nal measured data onto the basis functions (f>k(A) and reconstruct 
the unperturbed image using Eq. {JJ. One of the most interest- 
ing properties of this approach is that, since the real signal in 
each pixel only appears in the first few coefficients (a fundamen- 
tal consequence of the PCA decomposition), the influence of 
noise is largely minimized. W e deconvolve these image s using 
a Richardson-Lucy algorithm dRichardsonl 119721 : iLucvi 1 19741) 
controlling the number of iterations to avoid any ringing effect 
(although we hardly see any fringing given the absence of noise 
in the projection images). 



4. Inversion 



The SIR code (Ruiz Cob o & del Toro IniestaL[l992h was applied 
on the deconvolved spectra to derive the stratification of tem- 
perature T, magnetic field strength B, inclination y, azimuth <p, 
and line of sight velocity Vi os versus continuum optical depth r. 
The values of those parameters were retrieved at a number of 
optical depth points called nodes. We used a maximum of seven 
nodes in T, five in B, V\ os , y, and two in <p. Five different ini- 
tializations were used to avoid local minima. Because we are 
interested in the penumbra and our sunspot has negative polar- 
ity, we solved the 180° ambiguity of the azimuth by choosing the 
value in such a way that the magnetic field vector B points to the 
umbra. Additionally, we transformed B into the local reference 




1 Describ ed inlvan Noo rt (2012) and obtained from the Hinode pupil 
specified bv lSuematsu et al.l d2008t) . 
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Fig. 2. Filling factors of regions with reversed polarity at t — 0. 1 
after (solid line) and before (dashed line) deconvolution. 



frame. Likewise, we applied the same inversion procedure to the 
original data, i.e., before the deconvolution from the spatial PSF. 

5. Results and discussion 

In Fig. [TJwe present the reconstructed continuum intensity map 
of the region containing the sunspot, including the horizontal 
component of the magnetic filed vector in those pixels in which 
the field presents a reversed polarity. We note how the pixels 
showing reversed polarity are distributed over the whole penum- 
bra, with the exception of the most inner part. There is an obvi- 
ous concentration of reversed flux in some patches of the most 
external part of the outer penumbra. The full spot has a total 
unsigned magnetic flux of fl> z = 1.19 xlO 22 Mx, measured at 
t -0.1, while the flux of the penumbra is = 0.79 x 10 22 Mx. 
The reversed flux in the penumbra accounts for 8.4% of this 
penumbral flux, amounting to 17% of the pixels of the penum- 
bra, and 72% of these pixels present a positive velocity (down- 
flow). This figure should be compared with the 56% of the pix- 
els of the full penumbra harboring downflows. The correlation 
between downflows and reversed fields is probably tighter than 
this 72% because the velocities are measured in the line of sight 
reference frame, while B has been translated into the local ref- 
erence frame. The preceding figures increase to 30% (reversed 
flux over total unsigned flux), with 28% of the pixels harbor- 
ing reversed polarity and 88% of these presenting downflows at 
t =1. Nevertheless, B at r =1 has a large uncertainty (usually 
error bars around 600 G in B and 15° in y). 

In Fig.|2]we plot the ratio of pixels harboring inverse polarity 
with respect to the total number of pixels at t = 0.1 versus the 
radius from the center of the umbra, after and before deconvo- 
lution. The ratios shown in this figure are very similar to thos e 
found in numerical simulations (cf. Fig 13 in iRempeH 1201 2l) . 
except for the increase in the filling factor at radii larger than 
19 Mm. We need to take into account that our sunspot is not 
circular and that the upper left penumbra is certainly wider that 
the bottom right one. However, the similitude with the numeri- 
cal simulation results is clear, even considering the difference in 
spatial resolution. Also the inverse polarity pixels for the orig- 
inal data (before deconvolution) only appear at a radius larger 
than 14 Mm, while after deconvolution we already found inverse 
polarity pixels at 8 Mm. 

In Fig. [3] we show a zoom of two middle penumbra regions 
(marked with black rectangles in Fig. [TJ, one on the limb side 
(bottom panel) and the other on the center side (upper panel). 
It is conspicuous that the reversed field pixels are in dark fila- 
ments and, in many cases, they seem to mainly be placed at the 
boundary of bright filaments. The circles at (15,29) and (18,38) 
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Fig. 3. Sections in the center side (upper) and limb side (bottom) 
penumbra White arrows correspond to the horizontal compo- 
nent of the magnetic field in points with reversed polarity. White 
circles at (15,29) and (18,38) mark the position of the spectra 
shown in Fig. [4] 



mark the position of the pixels whose spectra are shown in Fig. 
|4] The first one is an example of a normal polarity pixel, and 
the second is a reversed one. The corresponding stratification of 
T, B, V z , and y in a local reference frame is plotted in panels 
in the bottom row. This pixel harbors an upflow and a negative 
polarity (y > 90°). The example of a reversed polarity pixel 
is plotted in the middle row. The corresponding stratification 
through the photosphere is plotted in the panels in the bottom 
row. The Stokes V profiles present a clear positive polarity (op- 
posite to that of the umbra). Also we note the small amplitude 
of these Stokes V profiles. The corresponding model is consider- 
ably cooler, and it presents a downflow and a weaker field with a 
inclination lower than 90° at all layers below optical depth 0.01. 
These profiles are contaminated by a small residual of the nor- 
mal polarity profile. There is a small reversed lobe at the blue 
wing of the Stokes V profiles of both 630. 1 nm and 630.2 nm 
Fe i lines. This small contamination and the consequent asym- 
metry of the Stokes profiles could be responsible for the strong 
gradient of both B and y, although the profiles are not perfectly 
reproduced by the model. 

It is well known that any overcorrection of stray light in a 
spectral image generates inverse polarity profiles at the border 
of bright structures. This is caused because a deconvolution is 
equivalent, in those pixels, to the substraction of a fraction of 
the spectral profile of the brighter neighbors. We can assess that 



our inverse profiles do not come from such an artificial effect 
for two reasons. First, we have undercorrected our data, because 
the used PSF is a lower limit of the real one (see, for instance, 
iJoshi et al.Ll2011h .As a consequence, the continuum contrast has 
only changed from 6.3% in the original data to 1 1.8% in the de- 
convolved data. This value is clearly below t he fig ure of 13-14% 
found in the simulations dDanilovic et al. L 120081) . Second, we 
checked that the original (nondeconvolved) Stokes V profiles of 
the 80% of those pixels that harbor an inverse polarity in the 
deconvolved map present the characteristic third lobe in the red 
wing of Stokes V in the original data (37% of the penumbr al pix- 
els ha rbor such a third lobe in the original data). Following Franz 
(201 lj), this is a hint of the presence of a hidden field with oppo- 
site polarity. With these two points we want to demonstrate that 
we are not introducing new information by the deconvolution of 
the PSF. We are simply disentangling mixed information. As an 
illustration, we present in Fig.|5]the Stokes profiles before the de- 
convolution of the same inverse polarity profile that we plotted 
in Fig. |4] Note the presence of the third, inverse polarity, lobe at 
the red wing of the Stokes V profiles of both spectral lines. We 
inverted these original (nondeconvolved) Stokes profiles using 
SIR, but now considering a two-component atmosphere. In the 
lower panels one can see the corresponding optical depth strati- 
fication of the different quantities for both models. Interestingly, 
the model of the red component is quite similar to the corre- 
sponding model of the deconvolved profile shown in Figure [4] 
The model of the blue component explains the contamination by 
stray light, and consequently it is similar to the model that re- 
produces a bright pixel. For this reason, the blue model is quite 
similar to the blue model presented in Figure[4] 

6. Conclusions 

We report here observations of reversed polarity fields at the 
border of bright penumbral filaments in the whole penumbra. 
The spatial distribution of the reversed polarity fields is similar 
to what is pred icted by the most recent numerical simulations 
dRempelL 120121) . This has been possible after deconvolution of 
the original data using a regularization method based on a prin- 
cipal component decomposition of the profiles. 

The Stokes profiles of those pixels that show a reversed po- 
larity in the deconvolved map display a small opposite polarity 
third lobe on the red wind in the orig inal data. Th is is a clear 
hint of hidden reversed polarity fields (Franzl l201 lh . We found 
a reversed flux in the penumbra at optical depth r = 0.1 of 8.4% 
of the total unsigned penumbral flux. It corresponds to 17% of 
the pixels of the penumbra. Seventy-two percent of these pix- 
els present positive velocity (downflow). All these facts indicate 
that we have found observational proof of convection in sunspot 
penumbra filaments. 

Acknowledgements. Hinode is a Japanese mission developed and launched by 
ISAS/JAXA, collaborating with NAOJ, NASA and STFC (UK). Scientific oper- 
ation of the Hinode mission is conducted by the Hinode science team organized 
at ISAS/JAXA. Support for the post-launch operation is provided by JAXA and 
NAOJ (Japan), STFC (U.K.), NASA, ESA, and NSC (Norway). Financial sup- 
port by the Spanish Ministry of Economy and Competitiveness through project 
AYA2010-18029 (Solar Magnetism and Astrophysical Spectropolarimetry) is 
gratefully acknowledged. AAR also acknowledges financial support through the 
Ramon y Cajal fellowship and the Consolider-Ingenio 2010 CSD2009-00038 
project. We thank C. Beck and M. van Noort for fruitful discussions. 

References 

Bellot Rubio, L. R., in Magnetic coupling between the Interior and the 
Atmosphere of the Sun, S.S. Hassan and Rutten (eds.), (Springer- Verlag: 
Berlin) ASP Ser, 2010, 193 



3 



B. Ruiz Cobo and A. Asensio Ramos: Returning magnetic flux in sunspot penumbrae 




8000 
7000 
6000 






-2 -1 

lOg T 




Fig. 4. Upper 4 panels: Stokes profiles and best fits obtained by the SIR inversion applied to a pixel in a penumbral filament 
harboring a magnetic field with the same polarity of the umbra. Middle 4 panels: the same for a pixel with the reversed polarity. 
Lower 4 panels: corresponding atmospheric models. The thick black line corresponds to the profiles shown in the upper row of this 
figure and the red ones to the middle row. 
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Fig. 5. Upper 4 panels: original (before deconvolution) Stokes profiles and best fits obtained by a 2-component SIR inversion applied 
to the inverse polarity pixel shown in Fig. [4] The two components are depicted in blue and red, and the combination in black line. 
Lower 4 panels: the two atmospheric models. 
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